Efficient Optimal Minimum Error Discrimination of 
Symmetric Quantum States 
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This paper deals with the quantum optimal discrimination among mixed quantum states enjoying 
geometrical uniform symmetry with respect to a reference density operator po. It is well-known that 
the minimal error probability is given by the positive operator-valued measure (POVM) obtained 
as a solution of a convex optimization problem, namely a set of operators satisfying geometrical 
symmetry, with respect to a reference operator IIo, and maximizing Tr(pollo). In this paper, by 
resolving the dual problem, we show that the same result is obtained by minimizing the trace of a 
semidefinite positive operator X commuting with the symmetry operator and such that X > po. 
The new formulation gives a deeper insight into the optimization problem and allows to obtain 
closed-form analytical solutions, as shown by a simple but not trivial explanatory example. Besides 
the theoretical interest, the result leads to semidefinite programming solutions of reduced complex- 
ity, allowing to extend the numerical performance evaluation to quantum communication systems 
modeled in Hilbert spaces of large dimension. 



I. INTRODUCTION 

In a quantum system Alice prepares the quantum chan- 
nel into one of several quantum states. Bob measures the 
quantum channel by a set of measurement operators and, 
on the basis of the result, it guesses the choice made by 
the transmitter. These actions lead to a classical chan- 
nel, and the problem arises of finding the measurement 
operators that provide optimal performance according to 
a predefined criterion, in this paper the minimum error 
probability. However, to solve the problem, except for 
some particular cases, has appeared to be a very diffi- 
cult task since the pioneering contributions in the seven- 
ties [H-i|. 

In recent years, particular attention has been paid to 
quantum states satisfying geometrical symmetry |^-u| ; m 
view of applications to optical communication systems. 
In some specific cases, including symmetric pure quan- 
tum states 0, Q and symmetric mixed quantum states 
with a characteristic structure Q, the solution named 
square root measurement (SRM) proves optimal. Never- 
theless, in general, SRM represents a suboptimal strat- 
egy, although it provides pretty good performance in 
many scenarios. 

In this paper, we are concerned with the construction 
of optimal POVM for the discrimination of symmetric 
mixed quantum states. We present some results that 
provide intuition into the problem and offer perspectives 
on its solution. 

Optimal quantum state discrimination represents a 
convex optimization problem, and, as such, it can be 
formulated in a primal and in a dual form, with the lat- 
ter having a reduced number of variables and constraints 
[H, [f| . In this paper we investigate how primal and dual 
problems simplify with symmetric quantum states. For 
the primal problem this study was already considered in 
Q. Herein, we extend the analysis to the dual prob- 
lem, where the optimal solution can be searched in a 



set of smaller dimension. The simplified formulation of 
the dual problem is illustrated with an example where 
a closed-form solution is easily found. For problems of 
large dimension, that cannot be solved analytically, the 
reduced number of variables in the simplified dual state- 
ment becomes useful to the numerical solution by means 
of semidefinite programming (SDP) tools. 



II. GENERAL FORMULATION OF MINIMUM 
ERROR DISCRIMINATION 

The quantum decision problem is formalized in a TV- 
dimensional complex Hilbert space JC [l|, where an en- 
semble of quantum states pi, i = 0, 1, . . . , M — 1, with 
prior probabilities qi > 0, Xa=o 1 <Zi 1 , is given. The 
quantum states pi are density operators on 5f, i.e., (self- 
adjoint) positive semidefinite (PSD) operators (p, > 0), 
with unit trace, Tr (pi) = 1. For notation convenience, we 
denote by T the class of the PSD operators on 3~C. The 
measurement operators n^, i = 0, 1, . . . , M— 1, constitute 
a POVM having the properties € CP and Yli=o * ^» = 
/ , where / is the identity operator on "K. The transition 
probabilities of the resulting quantum channel become 
p = Tr (pjlj) : so that the probability of correct de- 
tection is given by P c = J2i=o 1 (/°i n i) ■ 

Hence, the problem of finding the maximum proba- 
bility of correct state discrimination can be concisely 
stated as follows. 

Primal problem (PP1). Find the maximum of the 
probability of correct detection P c = J2i=o 1 9iTr (piUi) 
over the class of the POVM on JC. 

The analytical solution of PP1 is in general difficult 
since P c has to be maximized over the whole M-tuple of 
measure operators n^. As a matter of fact, closed- form 
results to the primal problem are available only for some 
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particular quantum mechanical systems, e.g., the binary 
case pj. Nevertheless, since the objective is to search a 
global maximum of a linear function into a convex set, 
the problem can be faced by means of numerical tools 
such as SDP. Besides, according to classical results in 
convex optimization theory in place of the primal 
problem, it is in general more convenient to consider its 
corresponding dual problem, since it presents a smaller 
number of variables and constraints [3|, [6J . 

Dual problem (DPI). Minimize the trace of the 
optimization operator X over the class CP, subject to the 
constraints X > QiPi, i = 0, 1, . . . , M — 1. Once found a 
minimum trace operator X opt , its trace gives the maxi- 
mum probability of correct detection, P c = Tr (^ pt)- 

The equalities (X opt - qtPi)^t = Hi(X opt - qipi) = 
0, i = 0, 1, . . . , M — 1, are necessary conditions on the 
optimal POVM. These conditions become sufficient, once 
the searched measure operators are constrained to belong 
to y and to solve the identity on "K. 



III. DISCRIMINATION OF SYMMETRIC 
QUANTUM STATES 

In quantum detection an important role is played by 
geometrically uniform symmetry Among the sev- 

eral generalizations, we consider the basic case of sym- 
metric mixed quantum states generated from a reference 
density operator po as 



S % po S~ 



0,1,...,M-1 



(1) 



where the symmetry operator S is unitary (SS^ = S^ S = 
I) and such that S M = I. The geometry implicitly 
requires that the mixed states are equiprobable, i.e., 
qi = 1/AI, i = 0, 1, . . . , M — 1. In Q it was shown that 
optimal POVM having the same symmetry can always be 
found. Hence, without loss of generality, we can assume 
that 



Ui = s i n s~ 



i = 0,l,...,M-l 



(2) 



where ITo is the reference measure operator. Conse- 
quently, for a fixed M, the knowledge of po, IIo and S 
is sufficient to fully describe the state ensemble and the 
POVM. 

The density operators (fTJ) have all the same rank, and 
the same holds for the measure operators ([2]). As proved 
in 0, the optimal measure operators can be assumed 
to have rank no higher than that of the corresponding 
density operators, namely rank (IIo) < rank(po) . 



A. Primal Problem for Symmetric Quantum States 

The specific geometry of the state ensemble can be 
exploited to get insight into how solving the state dis- 
crimination problem. In particular, the primal problem 
PP1 can be rewritten in a simpler form as follows. 



Primal problem for symmetric quantum states 
(PP2). Find the maximum of P c = Tr (p a Il ), with 

n 9 e? and such that J2fL 1 S l U S~ l = I. 
Proof: This formulation was first given in It 
can straightforwardly be proved by using (fTJ) and 
© in PP1, so that P c = E"o 1 * Tr (/'i n i) = 
EtLo 1 I? Tr (S'poS-* S l U S- 1 ) - Tr (p U ). Moreover, 

if n > o then n, = smoS-' 1 > o . ■ 



B. Dual Problem for Symmetric Quantum States 

The optimization problem PP2 is comparatively sim- 
ple, and, perhaps, this is the reason why no particular 
attention has been paid in the literature to the study of 
the dual theorem to obtain an alternative formulation. 
In the following we investigate this point. 

Dual problem for symmetric quantum states 
(DP2). Minimize the trace of the optimization operator 
X over the class 7, subject to the constraints X > j^po 
and XS = SX. Once found a minimum trace operator 
X op t, its trace gives the maximum probability of correct 
detection, P c = Tr (X opt ). 

Proof: Define p[ = S^-^p^S- 1 . Let V be the fea- 
sible set according to the general dual problem in Sec- 
tion[TTl i.e., the set of PSD operators X such that X > p' i7 
i = 0, 1, . . . , M - 1, and let V be the set of PSD opera- 
tors X' such that X' > p' and X'S = SX'. The proof is 
organized in two steps. In the first step it is shown that 
V C V, while in the second step it is proved that for any 

X e V there exist X' e V such that Tr(X') = Tr(X). 
Then, the search of X can be confined into V. 

- Step 1: If X' € V, for the commutativity between 

X I and S we get X' = SX'S- 1 and recursively X' = 
S % X'S~ l . Hence, for any i 



X' - p> =x>~ s* P ' s- 1 

= S^'S-' 1 - Stp'oS-* 

= s'ix' - p'^s- 1 > 

since X' > p' for assumption. 

- Step 2: For each X G V we consider 

M-l 



(3) 



x' = — V s~ l xs 

M f-f 



i=0 



Being X > p' i for each i, it follows that 



M-l 



i l 

- Po 



Moreover, recalling that S = I 

sx's- 1 = — y s-^-^xs 1 - 1 = X' 
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and then X' commutes with S. Finally, 
Af-l Af-l 
Tr(X') = jjJ2 Tr{S- l XS l ) = ^ Tr (X&S-*) 

i=0 i=0 
Af-l 

(=0 

and the proof is complete. I 

Therefore, the search of the unknown optimization op- 
erator X can be restricted to the subclass of 7 composed 
by the PSD operators that commute with the symmetry 
operator S. 

The optimal IIo can be found by the relations (X opt — 
1 g/9 )n = n (X op t - jjpo) = subject to IIo £ P and 

Y^ 1 s i u s- i = l. 

In the next section, we develop a formulation, where 
the commutation condition XS = SX is replaced by an 
alternative constraint. 



C. Alternative Formulation of the Dual Problem 

Since the symmetry operator S is known a-priori, 
it is possible to exploit its spectral characterization to 
rewrite the dual problem DP2 as follows. 

Other form for the dual problem for symmetric 
quantum states (DP3). Let Ao, Ai, . . . , Ajy_ x be the 
N < N distinct eigenvalues of S, and Ni be the multi- 
plicity of Ai . Let U be a basis of eigenvectors of S, with 
the first Nq eigenvectors corresponding to Ao, the next 
N\ eigenvectors corresponding to Ai, and so on. Then, 
minimize the trace of the optimization operator X over 
the subclass of 7 consisting of block- diagonal operators 
with blocks of dimension Ni, under the constraint 
X > jjW poll. Once found a minimum trace operator 
X opt , its trace gives the maximum probability of correct 
detection, P c = Tr (X opt ). 

Proof: By DP2 the optimization operator X can 
be chosen to commute with S. Therefore, given an 
eigenbasis U for S, we can write X = UXW where, by 
well-known results on simultaneous diagonalization of 
commuting self-adjoint operators [10], X turns out to be 
block diagonal with the size of the blocks given by the 
multiplicity of the eigenvalues of S. Moreover, we find 
that Tr (X) = Tr (UXW) = Tr (X) and the constraint 
X > jjpo becomes X > j^WpoU. ■ 

The commutative requirement XS = SX is then 
replaced by fixing the block-diagonal structure on X . 
Given an optimal X opt , the reference optimal operator 
ITo is solution of (UX opt W - jjPo^o = H (UX opt W - 

■jjpo) = 0, subject to n e 7 and J^fLo 1 s * n o S -< = I ■ 
The new formulation of the dual problem is given in a 
form that is particularly suitable for SDP computational 
tools [ll| and, moreover, it is analytically more tractable 
that the previous version. We now quantify the complex- 
ity of the different approaches. 



Table I: Number of real decision variables d, equality 
constraints C e and inequality constraints d for the opti- 
mization problems. 





d 


C e 


Ci 


PPI 


MN 2 


1 


M 


PP2 


N 2 


1 


1 


DPI 


N 2 





M 


DP 2 


N 2 





2 


DP 3 


T N N 2 





1 



D. Problem dimension and number of constraints 

The complexity of a linear program, to get the solution 
of a convex optimization problem, is hard to evaluate in 
terms of arithmetic operations (see Q for details). Nev- 
ertheless, since we are seeking a global optimum, the di- 
mension of the feasible region for the objective function 
gives an order of the complexity of the problem Q . 

The set of self-adjoint operators on Oi forms an 7V 2 - 
dimensional real vector space. Therefore, for a given op- 
timization problem we can find the dimension d of the 
correspondent real space on which the considered objec- 
tive function is defined. In other terms, we find the num- 
ber of real variables d in a given objective function. The 
results are summarized in Table U where C e and Ci repre- 
sent, respectively, the number of equality and inequality 
constraints for a given problem. Note that the PSD con- 
dition on self-adjoint operators is counted as an inequal- 
ity, e.g., the relation IIo £ CP is counted as the inequality 

rio > o. 

The dual problem DP3 presents the smaller number of 
variables and constraints among the considered optimiza- 
tion problems and, in particular, d — X^o* IS smaller 
than N 2 , depending on the spectrum of the symmetry op- 
erator S. If S has all distinct eigenvalues, then N = N, 
Ni = 1 for each i, and the optimization operator X in 
DP3 becomes diagonal, giving d~N< N 2 . 



IV. EXAMPLE OF APPLICATION 

In this section, the previous results are applied to the 
discrimination of an ensemble of M symmetric mixed 
quantum states on a 2-dimensional (TV = 2) complex 
Hilbert space. The symmetry operator is 



(4) 



and it represents a linear transformation given by a coun- 
terclockwise rotation through angle n/M. We assume 
that the reference density operator has the general form 



Po 



ft 



o 



(5) 
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where a and (3 are real numbers. Since po is PSD, the 
feasible values of a and j3 are constrained as < a < 
1 and |/3 1 < \/a(l — a), respectively. Without loss of 
generality we assume a > 1/2. 

The operator S has two non-degenerate eigenvalues 
equal to Ai(5) = e t7r / M and X 2 (S) = . There- 

fore, the corresponding eigenvectors define the basis 



U 



1 

71 



i 



i 



(6) 



Consequently, for the dual problem DP3 the optimization 
operator X has to be diagonal 



X 



■>-i 






$2 



(7) 



where, both x\ and x 2 are real and non-negative, being 
X PSD. The constraint X > jjW p U reads 



X > 



1 

2M 



1 

(2a- l)-i2/3 



(2a - 1) + i2/3 
1 



, (8) 



and after some simple algebra, 1 we find that it can be 
rewritten as 

(2Mxi - \)(2Mx 2 - 1) - [(2a - l) 2 + (2/3) 2 ] > , (9) 

with xi > l/(2M) and x 2 > l/(2M) . Hence, the min- 
imum of Tr (X) = £i + x-2 is obtained for x\ = x% = 
1/(2M)(1 + ^(2a - l) 2 + (2/3) 2 . In conclusion, the min- 
imum error probability P e = 1 — P c = 1 — Tr (X) is 



P e = — - — J(2a-\) 2 + (2(3f 



(10) 



We note that the first term on the right-hand side of 
(|10P corresponds to a blind guessing on the equiproba- 
ble elements belonging to the quantum state ensemble, 
while the second term represents the gain due to the op- 
timal quantum discrimination. It is interesting to ob- 
serve that the minimum error probability is obtained for 
/3 = y/ a(l — a) and it results 



Pr 



2 

l -Ti 



(11) 



For instance, (fTTj) holds for a — 1 and f3 = which is the 
case of study considered by Helstrom [l| and Ban et al. 

which models linearly dependent spin-1/2 quantum 
states, where the generating density operator has rank- 
one po = (V'oXV'o) with pure state (ipo) = ^ . When 
po is diagonal, i.e., f3 = 0, (|10|) simplifies as 



P e = 1 - a— 
M 



(12) 



The optimal reference measure operator IIo can 
be found from the conditions given in Section MI CI 

(ux opt W - ^po)n = n (ux opt W - = o, that 

in this example simplify as (x\ I — ■jjPq)YLq = T\q(x\ I — 
jjPo) = 0, being X opt = x\I. Note that these condi- 
tions also imply that polio = IloPo ■ The optimal IIo c & n 
then be numerically found by solving a linear system of 
equations including the additional requirements IIo S CP 

and £f£o 1 s ' 1 n o S~ l = I . 

It is useful to observe that the condition 
Y^fLo 1 S l Ho = I implicitly fixes the value of 
the trace of n . In fact, Tr (E^o 1 & n o s ~ l ) = 



E,=o lTr ( n o s- l s i ) 



E-io'Tr^no^) 

MTr (n ) and being Tr (I) =N it follows that 
Tr(n ) = N/M . By simple algebra it can be found a 
closed-form expression for IIo in the following two cases. 

1 . /3 = . The optimal IIo is given by 



2 

Jl 



1 




(13) 



In particular, setting a = 1/3 the same numerical 
results obtained in 11 211 are found. 



2. j3 = \/a(l — a) . The measure operator results 



IT 



The constraint 



Al 



(14) 



Tin S 1 = I becomes 



,o " Al 0' 

JT Ei=o Pi = I showing that the quantum state 
ensemble has a particular structure. Indeed, such 
a specificgeometry has been considered by Yuen 
et al. in [3|, (IV. 4)] and the results therein reported 
are in agreement with (fTTj) and (JTJJ) . 

The proposed formulation of the dual problem has also 
proved useful to numerically solve systems of large dimen- 
sions, where the computational complexity sets a severe 
limit to the possibility of finding an optimal solution. 
This is the case of pulse position modulation (PPM), a 
modulation format candidate for deep space communica- 
tions [lH- In quantum PPM the states are defined in a 
Hilbert space given by the tensorial product of M sub- 
spaces, each of dimension n, and, therefore, the overall 
space dimension N grows exponentially with the PPM 
order M, being N = n M 0. In [ij], DP3 is applied to 
quantum PPM using the software CVX for SDP [lj. The 
solution of DP3 results considerably faster than DPI. As 
a limit case, for values of TV about 1000 the discrimina- 
tion problem was successfully solved with DP3 while, on 
the same processing unit, CVX fails with DPI, because of 
computer memory limits. 



CONCLUSIONS 



It is recalled that an Hermitian matrix is PSD if and only if its 
principal minors are all non-negative. 



We have studied the dual problem for minimum error 
probability discrimination of symmetric quantum states. 
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It has been shown that the optimization operator, in the 
objective function, can be assumed to commute with the 
symmetry operator. This result leads to an alternative 
formulation of the dual problem, that presents a reduced 
number of variables and constraints. The obtained dual 
statement is convenient to find analytical solutions to 
the discrimination problem, as we showed with an illus- 
trative example. On the other hand, the new formulation 
also permits a computationally efficient numerical solu- 
tion by means of SDP methods. This property is partic- 
ularly useful to study the performance limits of quantum 
mechanical systems described by geometrically uniform 



quantum states on Hilbert spaces of large dimensions, 
such as modulated coherent states in optical communi- 
cations. 
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